Baryon bias and structure formation in an accelerating universe 
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In most models of dark energy the structure formation stops when the accelerated expansion 
begins. In contrast, we show that the coupling of dark energy to dark matter may induce the 
growth of perturbations even in the accelerated regime. In particular, we show that this occurs in 
the models proposed to solve the cosmic coincidence problem, in which the ratio of dark energy 
to dark matter is constant. Depending on the parameters, the growth may be much faster than 
in a standard matter-dominated era. Moreover, if the dark energy couples only to dark matter 
and not to baryons, as requested by the constraints imposed by local gravity measurements, the 
baryon fluctuations develop a constant, scale-independent, large-scale bias which is in principle 
directly observable. We find that a lower limit to the baryon bias b > 0.5 requires the total effective 
parameter of state We = 1 -\-p/p to be larger than 0.6 while a limit b > 0.73 would rule out the 
model. 
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I. INTRODUCTION 

The epoch of acceleration which the universe seems 
to be experiencing |Q is commonly regarded as a bar- 
ren ground for what concerns structure formation. In 
fact, during an accelerated expansion gravity is unable to 
win over the global expansion and the perturbations stop 
growing. Mathematically, this is seen immediately from 
the equation governing the evolution of the perturbations 
in the Newtonian approximation in a flat universe: 



1 



2 







(1) 



where H — dloga/dr is the Hubble constant in a confor- 
mally flat FRW metric ds"^ — a^{~dT'^ + Sijdx^dx^), the 
subscript c stands for cold dark matter (here we neglect 
the baryons) and the prime represents derivation with 
respect to a = log a . When the dark energy field re- 
sponsible of the acceleration becomes dominant, flc ^ 
and the dominant solution of Eq. (|^) becomes Sc ''^ const. 
Only if gravity can overcome the expansion the fiuctua- 
tions are able to grow. It appears then that to escape the 
sterility of the accelerated regime is necessary to prevent 
the vanishing of flc- 

As it has been shown in ref. , an epoch of accelera- 
tion with a non-vanishing Qc can be realized by coupHng 
dark matter to dark energy. In fact, a dark energy scalar 
field (f) governed by an exponential potential linearly cou- 
pled to dark matter yields, in a certain region of the pa- 
rameter space, an accelerated expansion with a constant 
ratio ric/^ip and a constant parameter of state w^, re- 
ferred to as a stationary accelerated era. Similar models 
have been discussed in [||, ^. In |^ we showed that in fact 
the conditions of constant fl^ and uniquely determine 
the potential and the coupling of the dark energy field. 
In this sense, the model we discuss below is the simplest 
stationary model: any other one must include at least 



another parameter to modulate the parameter of state. 
The main motivation to consider a stationary dynamics 
is that it would solve the cosmic coicidence problem 
of the near equivalence at the present of the dark en- 
ergy and dark matter densities 1^, 0| . The stationarity 
in fact ensures that the two components have an identi- 
cal scaling with time, at least from some time onward, 
regardless of the initial conditions. In |^ it was shown 
that, by a suitable modulation of the coupling, structure 
forms before the accelerated era. Further theoretical mo- 
tivations for coupled dark energy have been put forward 
in ref. |^. 

As it will be shown below, the coupling has three dis- 
tinct, but correlated, effects on Eq. (|^): first, as men- 
tioned, it gives a constant non-zero flc in the accelerated 
regime; second, adds to the "friction" (1 -I- H' /H)S'^ an 
extra term which, in general, may be either positive or 
negative; third, adds to the dynamical term —^ilcSc a 
negative contribution that enhances the gravity pull. 

The dark energy coupling is a new interaction that al- 
ways adds to gravity (see e.g. 0)- The coupling to 
the baryons is strongly constrained by the local gravity 
measurements, so that we assume for simplicity that the 
baryons are in fact not explicitely coupled to the dark 
energy as suggested in [0 and, in the context of dark 
energy, in |Q |l^(of course there remains the gravita- 
tional coupling). This species- dependent coupling breaks 
the equivalence principle, but in a way that is locally 
unobservable. However, we show that there is an effect 
which is observable on astrophysical scales and that may 
be employed to put a severe constraint on the model. In 
fact, the baryon perturbations grow in the linear regime 
with a constant, scale-independent, large-scale bias with 
respect to the dark matter perturbations, that is in prin- 
ciple observable. Although the bias can be either larger 
or smaller than unity, we find that all the accelerated 
models require b < 1 i.e. baryons less clustered than dark 
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matter (sometimes denotes anti-bias). Such a baryon bias 
would be a direct signature of an explicit dark matter- 
dark energy interaction, well distinguishable from most 
other hydrodynamical mechanisms of bias (see e.g. [^ ). 



II. COUPLED DARK ENERGY 

Consider three components, a scalar field (f), baryons 
and CDM described by the energy-momentum tensors 
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The CDM energy density parameter is obviously ilc — 
1 — — y2 — while we also have fi^ — x'^ + y'^, and 
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u^. The system is subject to the condition x'^ + 
y~ +u'^ < 1. 

The critical points of system (||) are listed in Tab. I. 
We denoted with We — I + Ptot/ Ptot = 1 + x'^ — the 
total parameter of state. On all critical points the scale 



T^ff.^): Tjiffb) and T^i,(c), respectively. General covariance factor expansion is given by a 
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i^, where 



requires the conservation of their sum, so that it is pos- 
sible to consider a coupling such that, for instance. 
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where = SnG, while the baryons are assumed uncou- 
pled, ^'^(f,).^ = because local gravity constraints indi- 
cate a baryon coupling (3^ < 0.01 H, 0. Let us derive the 
background equations in the fiat conformal FRW metric. 
The equations for this model have been already described 
in 1^ , in which a similar model (but with a variable cou- 
pling) was studied. Here we summarize their properties, 
restricting ourselves to the case in which radiation has 
already redshifted away. The conservation equations for 
the field 0, cold dark matter, and baryons, plus the Ein- 
stein equation, are 
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p = 2/(3we), while each component scales as a"^^' . In 
the table we also denoted g = 4/3^-1-4/3/1-1-18, and we used 
the subscripts b, c to denote the existence of baryons or 
matter, respectively, beside dark energy. In the same ta- 
ble we report the conditions of stability and acceleration 
of the critical points, denoting /i+ = {—f}+ ^/l8 + (P)/2. 

As it can be seen, the attractor a can be accelerated 
but flc 0, so that structure cannot grow, as in almost 
all models studied so far. Therefore, from now on we 
focus our attention on the global attractor be, the only 
critical point that may be stationary (i.e. ilc and fl^ 
finite and constant) and accelerated. On this attractor 
the two parameters /? and p, are uniquely fixed by the 
observed amount of Sic and by the present acceleration 
parameter (or equivalently by Wg — p/ip + P) )• For 
instance, flc — 0.20 and We — 0.23 gives p — 3 and 
/?= 10. 



III. DIFFERENTIAL GROWTH RATE 

Definining the perturbation variables S = 
Sp/p, ^ip — S(f>, ik^Sui/a — 9H, the following 
conservation equations for CDM, baryons and scalar 
field in the synchronous gauge for the wavenumber k are 
derived: 



where U{(j)) = [/ge~ v 2/3mk0 _ The coupHng /3 can be seen 
as the relative strength of the dark matter-dark energy 
interaction with respect to the gravitational force. The 
only parameters of our model are and p (the constant 
Uq can always be rescaled away by a redefinition of (j>). 
For /3 = /i = we reduce to the standard cosmological 
constant case, while for /3 = we recover the Ferreira & 
Joyce model |Q. As shown in ref. [Q, the coupling we 
assume here can be derived by a conformal transforma- 
tion of a Brans-Dicke model, which automatically leaves 
the radiation uncoupled. To decouple the baryons one 
needs to consider a two-metric Brans-Dicke Lagrangian 
as proposed in 

The system (g) is best studied in the new variables 

0, 111 X = 2/ = f yZT/S, andu= f Vp^TS 

. Then we obtain 
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% + W[-^^ + OcX^, (5) 

(6) 
(7) 



+ (2+:|-)</7'-f -|^(^+i/^'x + 2/i2y2<^ = /3r!c<^) 
Moreover we obtain the Einstein equation 

TTl 

h" = -{1 + —)h' - 2 {12lp'x ~ Qpy'^Lp) +i{5cVLc + 5bVLb). 
H 

(9) 



Now, deriving the 5'^ equation we obtain 



l + -|^-2/3x)<5^ + (^-l)|(5cr! 
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Table I: Critical points. 



Point X y u p We stability acceleration 
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-6(py'^H + (12 + 4/3^) yj'a; 

+ ^ip + ^h'x + ^i^y^ip). (10) 

In the small-scale Newtonian approximation we can take 
the limit k ^ oo. In Eq. this amounts to neglecting 
the derivatives of if and the potential term iJ.^y'^ip, which 
gives 

P 1 

— ^+-h'x^pnj,. (11) 

Substituting in Eq. ( |lo| ) and neglecting again ip' , (p" and 
the potential term, we obtain 

^" + + §^ - 2/3:e) - ^T^c^^c - ^Stilb = 0, (12) 

where 7 = 1 + 4/3^/3, and similarly for Si, 

^b + (^^ + jp)^'b-liSc^c + Sbnh) = o. (13) 

Eq. ( p!2| ) corrects the equation given in ref. |[l^, which 
had a wrong sign (the error gives only a minor effect for 
the small (3 considered in those papers). In Eq. ( p^ ) the 
differences with respect to Eq. (|^) that we mentioned in 
the introduction appear clearly: the friction term and the 
dynamical term are modified, and the value of Sic is con- 
stant along the stationary attractor. On the stationary 
attractor D,b —i- and Eqs. ( p^ and ( [l^ ) can be written 
as 

where x,We and flc are given in Table I as functions of 
the fundamental parameters fi, (3 for any critical point. 
The solutions are 5c = a™* and 5b = ba™'^ where 

m± = i [-4 + 3we + 4/3a: ± A] (14) 

b± = 3flc/(37r2c + 4/3xm±) (15) 

where = (247flc + (-4 -f 3u)e + 4/3a;)2) . The con- 
stant b = 5b /5c = 6+ is the bias factor of the grow- 
ing solution m = m+. The scalar field solution is 



ip « {H^a^P-^yP /kfSciPVLc + mbx). For small wave- 
lengths ip (which is proportional to 5p^/ ) is always 
much smaller than 5c, 5b at the present time (although it 
could outgrow the matter perturbations in the future if 
P>1). 

The solutions m±,b± apply to all the critical solu- 
tions of Table I. Let us now focus on the stationary 
attractor be- For /? = we recover the law m± = 
\ [-1 ± (2417c + 1)^^^] that holds in the uncoupled expo- 
nential case of Ferreira & Joyce |[l4| . Four crucial proper- 
ties of the solutions will be relevant for what follows: first, 
the perturbations grow (i.e. m > 0) for all the param- 
eters that make the stationary attractor stable; second, 
the baryons are anti-biased (i.e. 6 < 1) for the parame- 
ters that give acceleration; third, in the k ^ H limit (and 
in the linear regime) , the bias factor is scale independent 
and constant in time; and fourth, the bias is independent 
of the initial conditions. Numerical integrations of the 
full set of equations (@-^ that confirm and illustrate the 
dynamics are shown in Fig. 1. Notice that, in the future, 
the perturbations will reenter the horizon because of the 
acceleration, so that the subhorizon regime in which our 
solutions are valid will not hold indefinitely. 

The species-dependent coupling generates a biasing be- 
tween the baryon and the dark matter distributions. In 
contrast, the bias often discussed in literature concerns 
the distribution only of baryons clustered in luminous 
bodies jl^. To observe the baryon bias one should take 
into account also the baryons not collapsed in galaxies, 
e.g in Lyman-a clouds or in intracluster gas. A measure 
of the biasing of the total baryon distributions is possible 
in principle but is still largely undetermined. Even taking 
the extremely simplified approach that the baryon bias 
coincides with the galaxy bias we are confronted with the 
problem that the galaxy biasing depends on luminosity 
and type |^^. So we considered only very broad limits 
to b: since the acceleration requires anti-bias, we assume 
0.5 < 6 < 1. For a comparison, the likeHhood analysis of 
ref. gives b G (0.8, 1.9) (for IRAS galaxies, at 99% 
C.I.), but is restricted to ACDM models with primordial 
slope n — 1; other estimations do not exclude anti-bias, 
and might even require it Jl8| . 

In Fig. 2 we show all the various constraints. To 
summarize, they are: a) the present dark energy density 
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Figure 1: Numerical evolution of the density contrast for a 
100 Mpc/h perturbation of dark matter (continuous lines), 
baryons (dashed lines) and scalar field (dotted lines). Thick 
lines: f},/! — 1.5,3 (or fi^ = 0.55, We = 0.67), resulting in a 
bias b « 0.3. Thin lines: (3,ii = 0.25,3 (or fi^ = 0.5, We = 
0.92): here the dark matter and baryon curves are almost 
indistinguishable since 6 ~ 1. 



baryon bias, the trend is that more massive objects are 
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Figure 2: Constraints on the stationary model: below the 
horizontal line the expansion is accelerated; in the light grey 
region the bias is between 0.5 and 1; between the vertical lines 
ilcjy is within the observed range. The dark grey region is the 
surviving parameter space. 



0.6 < < 0.8; 6) the present acceleration {we < 2/3, 
implying (3 > /i/2); c) the baryon bias 1 > b > 0.5. On 
the stationary attractor there is a mapping between the 
fundamental parameters fi, (3 and the observables We, 
so one can plot the constraints on either pair of variables. 
It turns out that these conditions confine the parameters 
in the small dark shaded area, corresponding to 

We e (0.59 - 0.67), (16) 
or /3e (1.1 -1.4) , (2.0 -2.6). (17) 

Therefore, the parameters of the stationary attractor are 
determined to within 20% roughly. It is actually remark- 
able that an allowed region exists at all. The growth rate 
m is approximately 0.5 in this region. For h > 0.73 the 
possibility of a stationary accelerated attractor able to 
solve the coincidence problem would be ruled out. If one 
considers the tighter limit < 0.4 for the supernovae 
la given at two sigma in ref. |^ for stationary attractors 
the allowed region would be further reduced, possibly re- 
quiring a lower b to survive. 

IV. CONCLUSIONS 

We have shown that if the universe is experiencing a 
stationary epoch capable of solving the cosmic coinci- 
dence problem then two novel features arise in the stan- 
dard picture of structure formation. First, a non-zero flc 
during the accelerated regime allows structure to grow; 
second, since the baryons have to be uncoupled (or very 
weakly coupled), the growth is species-dependent, result- 
ing in a constant baryon bias independent of initial con- 
ditions. Although there are no direct observations of the 



more biased with respect to the dark matter distribu- 
tion, so probably the galaxy bias is higher than the to- 
tal baryon bias. If this is correct, then b can be smaller 
than unity, as we find to occur for all accelerated models. 
We find that the bias strongly constrains the existence 
of a stationary epoch. Putting 6 > 0.5, and requiring 
0.6 < < 0.8, we get that the two free parameters /i 
and l3 are fixed to a precision of 20% roughly, while the 
effective parameter of state We is larger than 0.59. A 
higher bias or a lower We can easily result in ruling out 
this class of stationary models. On the other hand, the 
observation of a constant, scale- independent, large-scale 
anti-bias would constitute a strong indication in favor of 
a dark matter-dark energy coupling. 

The growth rate m is another observable quantity that 
can be employed to test the stationarity, for instance es- 
timating the evolution of clustering with redshift. So far 
the uncertainties of this method are overwhelming (see 
e.g. ||l^) but future data should dramatically improve 
its validity. The combined test of b and m will be a 
very powerful test for the dark matter-dark energy in- 
teraction. There is also the possibility to compare the 
resulting power spectrum or cluster abundance with ob- 
servations, although then we should know exactly when 
the stationarity begun, and what dynamics preceded it 
(see e.g. the variable coupling model of ref. |^). 

Although we investigated only the simplest stationary 
model, in which We is constant (a reasonable assumption 
over a small redshift range) , it is obvious to expect that a 
similar baryon bias develops whenever there is a species- 
dependent coupling; this, in turn, is requested to provide 
stationarity without conflicting with local gravity exper- 
iments. Therefore, we conjecture that the baryon bias is 
a strong test for all stationary dynamics. 
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